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We measure how strong, localized contact adhesion between grains affects the maximum static 
critical angle, C , of a dry sand pile. By mixing dimer grains, each consisting of two spheres that 
have been rigidly bonded together, with simple spherical monomer grains, we create sandpiles that 
contain strong localized adhesion between a given particle and at most one of its neighbors. We 
find that tan# c increases from 0.45 to 1.1 and the grain packing fraction, decreases from 0.58 
to 0.52 as we increase the relative number fraction of dimer particles in the pile, ua, from to 1. 
We attribute the increase in ta,n8 c (vd) to the enhanced stability of dimers on the surface, which 
reduces the density of monomers that need to be accomodated in the most stable surface traps. A 
full characterization and geometrical stability analysis of surface traps provides a good quantitative 
agreement between experiment and theory over a wide range of Ud, without any fitting parameters. 
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I. INTRODUCTION 

The presence of adhesive forces between grains can 
greatly alter the physical behavior of sandpiles jl],|| . Al- 
though the importance of intergrain adhesion has been 
noted in many fields, ranging from soil science to civil 
engineering, the understanding of the physical principles 
governing the link between the macroscopic behavior of 
sandpiles in which adhesion is present and the micro- 
scopic attractive force distribution within the sandpile 
remains limited. 

Some light has been shed on this subject recently 
through experimental and theoretical studies which have 
shown that small quantities of liquid added to a sandpile 
comprised of rough spherical grains can cause sufficient 
intergrain adhesion so that the angle of repose after fail- 
ure Q and also the maximum static angle of stability of 
the sandpile before failure Q, known as the critical angle, 
9 C , greatly increase. A continuum theory that links stress 
criteria for the macroscopic failure of the wet pile to the 
cohesion between grains, which in turn was attributed to 
the formation of liquid menisci with radii of curvature 
that were determined by the surface roughness charac- 
teristics of individual grains, has provided a satisfying 
quantitative explanation of the increase of C with the 
liquid volume fraction and air-liquid interfacial tension 




In this theory, each intergrain meniscus is assumed to 
exert the same average attractive force everywhere in the 
pile. This assumption appears to be valid for sufficiently 
wet sandpiles; however, when the volume of the wetting 
fluid is small enough this theory is incapable of explain- 
ing the data. There are two principal features of the 
small-fluid-volume data that are incompatible with the 
continuum theory ||. The first is that the increase in 
9 C with small amounts of wetting fluid is independent 
of the surface tension of that fluid [[||(| . The second is 
that, in order to quantitatively fit the data, one must 



assume that a small fraction of the wetting fluid is se- 
questered on the grains in such a way that it does not 
participate in the formation of inter-grain menisci that 
contribute to the adhesive stresses within the pile. It 
is possible that, at vanishingly small fluid coverage, the 
physical/chemical inhomogeneities of the grains' surface 
prevents the transport of the wetting fluid from one men- 
sicus to another thereby allowing a wide distribution of 
inter-grain cohesive forces. To begin to understand the 
effect of such a broad distribution of inter-grain forces 
upon the macroscopic properties of the sandpile, we con- 
sider in this paper an extreme example of such a distri- 
bution in which some inter-grain contacts have an arbi- 
trarily large cohesion while others have no cohesion at 
all. Consistent with the notion that the non-uniform dis- 
tribution of inter-grain cohesion is primarily significant 
at low fluid volumes, we study a system designed so that 
each grain has at most one strong cohesive contract. 

In order to begin to investigate how nonuniform dis- 
tributions of microscopic attractive forces between grains 
can affect the overall macroscopic stability of sandpiles, 
we have measured maximum stability angles of dry sand- 
piles made by mixing a weight fraction v& of dimer grains 
(two spherical grains rigidly bonded together) into spher- 
ical monomer grains. The measured \,m\Q c (vd) gradually 
increases over the entire range of < Vd < 1, despite a 
moderate drop in the total packing fraction of grains, 
within the pile, due to the more inefficient packing of the 
dimers. 

A detailed theoretical study of the failure mechanism 
of such piles leads us to the following key observations 
and conclusions, which enable us to quantitatively ac- 
count for the increase in tan0 c (z^): 

(i) For piles consisting of perfectly rough particles 
(large intergrain friction), the stability of the free sur- 
face upon tilting is limited by the particles on the surface 
layer, which fail by rolling out of the surface traps they 
sit in. 
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(ii) For a given surface trap geometry, dimer particles 
typically remain stable up to larger tilt angles; thus, for a 
mixture of monomers and dimers, pile stability is limited 
by the monomers on the surface, provided that the dimer 
concentration is not too large. 

(iii) Provided that individual grains rolling out of un- 
stable traps do not initiate avalanches, the pile will re- 
main stable as long as the density of stable surface traps 
is larger than the density of monomers on the surface 
layer. 

(iv) The ratio of the density of monomers on the sur- 
face layer to the total density of surface traps is {l—Vd)/2. 

(v) A statistical characterization of the particle-scale 
roughness of the surface associated with grain packing is 
necessary to determine 9 C quantitatively. 

The stability criteria for a pile consisting of a mix- 
ture of monomers and dimers with ideally rough surfaces 
can be cast as a purely geometrical problem under condi- 
tions where rolling-initiated surface failure is the primary 
mechanism that limits the stability of the pile. In this 
picture, monomers and dimers on the surface layer oc- 
cupy surface traps formed by the particles underneath, 
each of which have a different stability criterion associ- 
ated with the trap's shape and orientation with respect 
to the average surface normal and the downhill direction. 
For a pile to be stable at a given tilt angle, all the grains 
on the surface have to be sitting in a stable surface trap, 
suggesting that the least stable surface traps would con- 
trol the overall stability of the pile. 

For a pile consisting of monomers, there are actually 
twice as many surface traps as surface grains, and upon 
very gradually increasing the average tilt of the pile with- 
out disturbing the underlying grains, one finds that sur- 
face grains in traps that become unstable upon tilting can 
briefly roll down the pile's surface until they encounter 
an unoccupied stable surface trap which ends their de- 
scent, provided that at least half the surface traps remain 
stable so that all surface grains can be accomodated. 

A detailed analysis in Sec. [V shows that it is much 
easier to trap dimers than monomers on a perfectly or- 
dered close packed surface of spheres. This suggests that 
the stability of a random mixed pile of monomer and 
dimers is actually limited by the monomers on the sur- 
face, which accumulate in the most stable surface traps as 
9 C is approached. If one assumes that the dimers remain 
essentially stable, one might expect that only a fraction 
(1 — Vd)j e l of the surface traps remain stable when the 
pile is tilted to its maximum stability angle 9 c {vd)] at 
this angle, there are just enough to accomodate all the 
monomers on the surface layer. A detailed experimental 
characterization of the positions of s urfac e grains of ran- 
dom monomer and dimer piles in Sec. V B , combined with 
the computed stability criteria for m onomers and dimers 
occupying surface traps in Sec. VA, provide a quantita- 
tive explanation of the measured increase in 9 c (vd) with 
no adjustable parameters up to Vd ~ 0.6, where the as- 
sumption of monomer failure begins to break down. 

This good agreement over a wide range of Vd < 0.6 re- 



flects the subtle interplay between the distribution of ori- 
cntational and size fluctuations of the surface traps, some 
of which are stabilizing and some destabilizing. Thus, 
accurate characterization of "grain-scale" roughness (as 
distinct from the microscopic surface roughness of the 
grains) is essential to achieve quantitative agreement be- 
tween theory and experiment. The important role played 
by grain-scale roughness is also evident in the rheology 
of gravity-driven chute flows |?],|8| , where the precise na- 
ture of the bottom surface has significant influence on 
the resulting flow. Nevertheless, the results convincingly 
demonstrate that the stability of cohesionless grains with 
large inter-grain friction is indeed controlled by surface 
failure, and should otherwise be insensitive to the type 
of grain material. 

In addition to clarifying the role of surface failure in 
the stability of piles by connecting macroscopic measure- 
ments of stability angle to grain-scale composition, these 
results for a well-characterized sandpile also provide a 
critical link between the attractive forces within the sand- 
pile and the nonspherical geometry of a well-known frac- 
tion of constituent grains. In this respect, these mea- 
surements provide quantitative insight into similar mea- 
surements of the angle of repose after dynamic failure 
of less well-controlled piles of spheres and cylinders (e.g. 
peas and rice). From this perspective, dimers may be 
imagined as short elongated grains that have surface ir- 
regularities of the same order as the grain size. These 
irregularities promote the strong interlocking of adjacent 
grains, which inhibits the failure of the pile more than the 
typical contacts between smooth cylinders and ellipsoids. 

The rest of the manuscript is arranged in the follow- 
ing way. In Sec. we present the experimental method 
for preparing the dimer and monomer grains, measur- 
ing the critical angle of stability and grain packing frac- 
tion of the pile as a function of the dimer content, and 
characterizing the surface configuration of monomer and 
dimer piles. The results of these measurements are also 
reported in this section. The development of the theo- 
retical understa ndin g of these results starts in the next 
section, Section III, in which the stability of a pile [of 
spherical grains] is posed as a geometrical problem. This 
section summarizes an earlier attempt to treat pile sta- 
bility geometrically || and presents a different and more 
general strategy for the solution, which can be extended 
to include dimer grains. Using this approach, in Sec. |TV| , 
we fully solve the stability problem on a triangular close- 
packed surface layer for monomers and dimers. Section |v| 
presents a statistical analysis of the measured shape and 
orientation of traps on real surfaces of piles comprised 
solely of either mon omers or dimers, and uses the method 
outlined in Sec. Ill to determine the corresponding solu- 
tion to the stability angle. In Sec. VI, we summarize the 



main findings and insight gained from this study, as well 
as possible future directions. 
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II. EXPERIMENTAL 



completely of dimers. 



We prepare the dimer grains by bonding glass spheres 
of radius d = 4.6 ± 0.2 mm and density p g — 2.3 g/cm 3 
together using methyl acrylate glue. The glue in its car- 
rier solvent completely coats the surfaces of the glass 
spheres and accumulates in a contact meniscus between 
the spheres. After less than one day of drying, a strong 
shear-rigid bond between the spheres is formed. The vol- 
ume of the dried glue is much smaller than the volume of 
the spheres, so the two bonded grains have the appear- 
ance of an ideal dimer or doublet. Because the glue coats 
the entire surfaces of the grains and may thus alter the 
friction coefficient of the grains, we have likewise coated 
all the monomer grains with methyl acrylate so that the 
friction coefficient at the contact points between grains, 
whether monomers or dimers, is identical. The sandpiles 
are prepared by mixing together varying weight fractions 
Vd of dimer grains into monomers, as determined by a 
balance. The mixtures are placed in a clear plastic box 
having a square bottom that is 6.5 cm wide and a height 
of 5.5 cm, yielding an average number of spheres per pile 
of about three hundred. 

To measure the critical angle, we employ a procedure 
that is identical to one that was used to study the critical 
angles of wet sandpiles We tilt the box at an angle, 
as shown in Fig. |l](a), and shake it back and forth about 
five times along the direction of the lower edge (normal 
to the page in the figure). This distributes the grains so 
that the surface of the pile is normal to the direction of 
gravity, as shown in Fig. 0(b). There is no noticeable 
size segregation of the grains introduced by the shaking. 
We have purposefully avoided tapping the pile in order 
to prevent a densification that could affect the pile's sur- 
face characteristics and stability We then place the 
lower edge of the box on a table and slowly tilt it so 
that the bottom rests flat on the table and normal to 
the direction of gravity. If the pile fails catastrophically, 
then no measurement is recorded, but if the pile remains 
stable after only isolated movement and resettling of a 
few grains, then a measurement of the static angle of the 
pile, 0, is recorded, as shown in Fig. |l](c). After several 
trials, a rough determination of the stable angle is ob- 
tained, and thereafter the initial angles are not chosen 
randomly, but instead are kept close to this value. Us- 
ing the results of ten trials, we average the values of the 
three largest angles to obtain the critical angle, 9 C . This 
value of the critical angle is reproducible; the variation 
in the three angles used to obtain the average is about 
ten percent for all values of v&. This procedure yields 
slightly larger angles than those found in typical angle of 
repose measurements in which the sandpile is induced to 
fail. 

Figure || depicts 9 C as a function of v&. We find that 
it increases approximately linearly, from tan# c ss 0.45, 
a well established value for a wide variety of dry spher- 
ical grains 0, to tan# c ss 1.1 for a sandpile comprised 



We have also measured the average packing volume 
fraction of grains, <E>, in the sandpile as a function of v& 
by measuring the mass of water required to fill the voids 
in the pile as it stands in the tilted configuration shown 
in Fig. |l|(a). Based on the relative number of grains at 
the surfaces compared to those within the pile, we esti- 
mate that the lowering of the local grain packing fraction 
due to the open or wall surfaces makes the measured $ 
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FIG. 1. Measurement of the critical angle of stability, 
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FIG. 2. The critical angle of stability, 6 C , as a function of 
dimer weight fraction, Vd- 




FIG. 3. The packing fraction, <£>, as a function of dimer 
weight fraction, Vd- 
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appear to be about two percent smaller than the corre- 
sponding bulk value. These measurements of <&{vd) are 
plotted in Fig. |[ The overall reduction in the packing 
fraction of about ten percent indicates that dimer grains 
pack less efficiently than monomer grains. 

In order to experimentally characterize the surfaces of 
random piles of either monomer or dimer grains, we have 
taken stereo digital images of stable piles using top and 
front views in order to reconstruct the three-dimensional 
coordinates of the centers of all of the spheres visible on 
the surfaces of the piles. We do not include spheres that 
touch walls or the bottom surface of the container. Due 
to the significant surface roughness of the piles, especially 
in the case of the dimer pile, we occasionally detect the 
position of a sphere that lies more than one diameter be- 



low the average surface defined by all the grains. These 
spheres are not true surface spheres and are eliminated 
from conside ratio n in the subsequent surface trap analy- 
sis (See Sec. 



VB 



The measured positions of the surface spheres for a sta- 
ble monomer pile and a dimer pile are shown in Fig. |J(a) 
and |](b), respectively. The monomer pile shown is very 
close to its critical angle of stability (9 — 23°), whereas 
the dimer pile shown here, while at a much higher angle 
(6 = 41°), is still somewhat below its 6 C . Finally, we have 
qualitatively observed that the roughness of the sand- 
pile's free surface increases somewhat as more dimers are 
included in the pile. This increase has been quantified in 
terms of larger fluctuations in the shapes and o rientations 
of the surface traps, as presented in Sec. VB. 




FIG. 4. Reconstructed positions of spheres at the surface of a pile consisting of (a) monomers, and (b) dimers. The technique 
used can identify the spheres, but not the bonds of the dimers, thu s onl y half of a dimer might be shown in some cases. However, 
that information is not needed in the subsequent analysis in Sec. |VB|. 



III. PILE STABILITY AS A GEOMETRICAL 
PROBLEM 



A. Background and Context 

Various measurements of the angle of repose 9 r for 
cohesionless piles of smooth spherical particles come up 
with the same value of about 22°, largely independent of 
the makeup of the spheres or of their surface properties 
P|. On the other hand, the shape of particles in a pile 
has a large influence on 8 r and the slightly larger crit- 
ical angle 8 C |11|. Furthermore, our granular dynamics 
simulations of piles, made of spheres with Hertzian con- 
tacts and static friction, show that 9 r initially increases 
rapidly with increasing friction coefficient /i, although it 
saturates at a value of about 22° for fx > 1 (see Fig. ||). 

These observations suggest a primarily geometrical ori- 
gin for the robustness of 9 C for perfectly rough spheres 
(/i ^> I), which can be further studied in an idealized sys- 
tem in which sliding is disallowed, due cither to a very 



large friction coefficient or to interlocking surface irregu- 
larities. The spheres in such a static pile can be classified 
into two groups as follows: A surface layer that consists of 
spheres held in place by exactly three spheres and their 
own weight, and interior spheres that have more than 
three contacts. (We do not consider the small popula- 
tion of "rattler" spheres with three contacts that might 
exist in the interior of the pile, which will presumably 
not influence the stability of the pile.) The centers of 
the three spheres that support each sphere in the surface 
layer form the vertices of what we henceforth call the 
"base triangle" associated with that sphere. 

As the tilt angle of the pile is increased, spheres on the 
surface layer can move by rolling out of the surface trap 
formed by the three supporting neighbors. However, in- 
terior particles (excluding rattlers) are held in place by a 
cage formed by their contacting neighbors, and if the fric- 
tion coefficient is sufficiently large to preclude any sliding, 
they cannot move until this cage is destroyed by the mo- 
tion of at least one of their neighbors. This suggests that 
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initiation of failure occurs at the surface layer, provided 
that sliding is disallowed. If the coefficient for rolling 
friction is small enough to be neglected, the stability of a 
sphere on the surface layer, and consequently the deter- 
mination of 9 C , becomes a purely geometrical problem. 
For finite values of fj,, other failure mechanisms can be 
expected to reduce 9 C from this surface controlled value, 
as observed in Fig. [|. 

A recent attempt at a theoretical determination of 9 C 
from the perspective of surface stability was made by 
Albert and co-workers |jj . They have considered the sta- 
bility of spheres at the surface, supported by three close- 
packed spheres that form a base triangle, and calculated 
the tilt angle 9 ma x at which the sphere would roll out 
of the trap formed by the base triangle as a function of 
"yaw" <f>, i.e., the relative angle of orientation of the tri- 
angle with respect to the downslope direction: 



t&n9 max (4>) 



2v / 2cos(</>) 

9max(4>) = 9max{4> + 27r/3) 



101 < 



3' 



(1) 
(2) 



This stability criterion is periodic with period 27r/3 due 
to symmetry. Yaw </> = corresponds to an orientation in 
which one of the edges of the base triangle is perpendic- 
ular to the downslope direction. In order to account for 
the randomness in the orientations of the base triangles 
on a disordered surface, they have suggested that the ap- 
propriate value for 9 C can be obtained by averaging 9 max 
over yaw <fi; they assumed a uniform distribution for this 
quantity. This yielded a value for 9 C that closely matched 
experimental observations. 
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FIG. 5. Angle of repose for a pile of spheres with Hertzian 
contact interactions and static friction as a function of friction 
coefficient fi betweeen the spheres, obtained from granular dy- 
namics simulations by extrapolating flow rates as a function 
of tilt angle to zero flow. Details of the simulation technique 
can be found in Refs. MM. 



B. Approach 

As already pointed out by its authors, the calculation 
in Ref. |)| represents a mean-field approximation, since 
it ignores variations in the shape and orientation of indi- 
vidual surface traps, which can be parametrized by their 
local tilt, yaw and roll angles, and the actual edge lengths 
of the base triangles. (For definitions of the these param- 
eters, see the Appendix.) Nevertheless, their result is in 
good agreement with experiments. We will address the 
dependence of the stability of a surface trap on so me o f 



these additional parameters in more detail in Section V A 



There is, however, a more serious complication with 
the averaging approach than the neglect of fluctuation 
effects: The stability of the pile requires all particles on 
the surface to be stable. Thus, the stability of the pile 
should be dictated by the particle in the least stable sur- 
face trap, and not an average stability criterion. One 
might thus wonder why the averaging approach appears 
to work so well. 

In fact, for a pile of monomers, the number of base tri- 
angles forming potential surface traps is essentially twice 
the number of surface particles that actually reside in 
them. This is easy to see in the case of close-packed lay- 
ers, where each successive layer to be placed on top has a 
choice among two sublattice positions; this is what gives 
rise to random stacking. 

The relationship is actually more general. If the sur- 
face of a pile is sufficiently smooth such that an average 
surface normal vector to the pile can be determined (note 
that this is a prerequisite to actually being able to define 
and measure 9 C ), the base triangles associated with sur- 
face traps can be identified by a Delaunay triangulation 
of the sphere centers at the surface layer, projected onto 
the plane of the mean pile surface. In such a triangu- 
lation, the number of triangles per surface layer sphere 
is exactly two, since the sum of all the interior angles 
of the triangles is 7rx(no. triangles)= 27rx(no. surface 
particles). An intrinsic assumption here is that the sur- 
face layer is similar to the "sublayer" , consisting of those 
spheres that would become part of the new surface layer 
if all the original surface particles were removed simulta- 
neously. The triangulation procedure to identify the sur- 
face normal vector and all of the potential surface traps 



is discussed in greater detail in Sec. VB 



This ratio of surface trap to surface sphere density in- 
dicates that in a stable pile, only half of the traps are 
actually filled. The pile will then find a stable configura- 
tion as long as at least half of the surface traps are stable 
at the given tilt angle of the pile, since surface spheres 
that are in unfavorable traps can roll down the slope un- 
til they find a vacant trap of sufficient stability, assum- 
ing that they do not gain enough kinetic energy to knock 
other particles off their traps and cause an avalanche. 
Continuous failure of surface spheres will occur if there 
are never enough traps to stabilize the entire layer. 

This leads to the conclusion that the stability of the 
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pile is actually determined by the median stability an- 
gle of the traps, not the mean. Nevertheless, as shown in 
Sec. IV, the quantitative difference between this criterion 
and that studied in Ref. |J is small; about 1.6 degrees. 



IV. MONOMER-DIMER STABILITY ON A FLAT 
CLOSE-PACKED SURFACE 

Before launching a full-scale analysis of the pile sta- 
bility problem for random piles having random surface 
grain configurations, it is instructive to consider the im- 
plications and power of this new approach to stability on 
a simplified system consisting of a mixture of spheres 
(monomers) and dimers sitting on a triangular close- 
packed lattice. The stability analysis leading to Eq.(Q) 
can be extended to dimers. Dimers sit in surface traps 
such that the vector connecting the centers of the two 
spheres forming the dimer are always parallel to one of 
the edges of the base triangles, thus their orientation with 
respect to the downhill direction can be described by the 
angle cj) as well. The resulting stability angle as a function 
of (j>, defined in the interval (— it, it), is: 



tanf 



r (0) 



2\/2cos(» ' 
V2 

cos(0-27r/3) : 

1 

V2cos(0-7r) ' 



\(f>\ < arctan(3v / 3), 
arctan(3\/3) <\<j>\< 2tt/3, 
2tt/3 < |0| < 7T. 

(3) 



The two functions defined by Eqs.(Q) and (||) are plot- 
ted in Fig. H(a). It is striking how much more stable 
dimers are compared to spheres for certain orientations 
of the traps. This is due to the more favorable position 
of the center of mass of the dimer, located where the two 
spheres meet, which makes it more difficult to roll out of 
the surface traps. An alternate, and perhaps more vivid, 
way of seeing the relative stability of dimers with respect 
to monomers is to plot the fraction of stable traps / s tab(#) 
at a given tilt angle 8: 



,77/2 

/ stab (0)= / de'Dos{e'), 

Je 



DOS (8)= / d<f>P ya , w (4>)5 {6 ■ 

J — 7T 

Vy aw ((p) = — , -7T < (j) < I"- 



2tt 



(4) 

(5) 

(6) 
(7) 



In the above, / s tab(#) is defined in terms of the den- 
sity of surface traps at a given stability angle, DOS (9). 
For this particular case, the distribution of yaw 7- > yaw ((/>) 
is assumed to be uniform in the interval (— tt, 7t), corre- 
sponding to an isotropic surface geometry. The result- 
ing plot of /stab(^) for monomers and dimers is shown in 
Fig. ||(b), and clearly demonstrates the difference in their 



stability. Consequently the critical angles of stability in- 
ferred from the median stability angle for monomers and 
dimers are: 



□monomer 



= 22.2° 

jdimor = 3^30 



(8) 
(9) 



Values obtained through the averaging procedure of Al- 
bert et d. |§ are only slightly larger: 6>™ onomcr = 23.8° 
111 and 8f mcl = 39.5°. 
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FIG. 6. (a) Stability angle of monomers (dashed line) and 
dimers (solid line) on a close-packed surface as a function of 
orientation (yaw) <f). Monomer and dimer curves are identi- 
cal for < (f> < 7r/3. (b) The fraction of stable traps at a 
given tilt angle 8 corresponding to a population of monomers 
(dashed line) and dimers (solid line) in surface traps with a 
uniform yaw distribution. 

For a surface layer consisting of a mixture of spheres 
and dimers, with a dimer volume fraction of v^, the sta- 
bility of the surface layer will be primarily controlled by 
the monomers, since dimers will be stable at most loca- 
tions on the surface at C and do not need to be consid- 
ered as surface particles for the purposes of the stability 
analysis. Thus, for a given tilt angle, the surface layer 
can find a stable configuration as long as the fraction of 
stable traps at that angle, / s tab(0), exceeds (1 — i^)/2, 
the density needed to accomodate all the monomers. 

The sequential filling of surface traps starting from the 
most stable one is somewhat analogous to the filling of 
an energy band in a fermionic system, with —9 ma x for a 
trap corresponding to the energy E of a fermionic state. 
DOS{9) can then be interpreted as a "Density of States" . 
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The monomer pile is analogous to a half-filled energy 
band and the addition of dimers lowers the filling frac- 
tion from 1/2. Thus, the critical stability angle 9™ lx (i/d) 
is determined by the "Fermi energy" of the system at the 
given filling fraction, defined through the implicit relation 

1 - v d 



/stab (C iX M) = 



(10) 



This relation has been plotted as a dashed line along with 
experimental data for 9 C in Fig. fj]. Although it captures 
the essential features of the dependence on dimer mass 
fraction and provides a compelling mechanism for this ef- 
fect, the results are not quantitatively comparable. The 
origin of the discrepancy lies primarily in the simplifica- 
tions made in characterizing the surface: Fluctuations in 
the shape and orientation of surface traps will broaden 
the DOS spectrum and consequently change the values 
obtained for 9 C . In Section [V|, we include the most rel- 
evant of such fluctuations in the analysis and compare 
results to available experimental data on the properties 
of such surfaces. Agreement between theory and experi- 
ment improves significantly when such fluctuation effects 
are taken into account. 
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FIG. 7. Measured and calculated values of 6 C as a func- 
tion of dimer weight fraction Vd- Circles: Experiment. 
Dashed Line: Computation based on a randomly oriented, 
flat, close-packed surface. 



V. MONOMER STABILITY ON A RANDOM 
SURFACE 

In generalizing the approach of the previous section to 
random surfaces, we will assume that the stability of the 
pile is still controlled entirely by the monomers. This 
assumption is likely to break down at very large dimer 
concentrations, so that comparison to experiments may 
not be appropriate in that case. However, it enables the 
determination of stability angles based entirely on the 
behavior of monomers, and avoids the extremely tedious 
analysis of dimer surface traps. Monomer traps, on the 
other hand, are completely characterized by their base 
triangle, formed by the centers of the three supporting 
spheres. 

There are two steps that are needed to obtain the 
monomer DOS required to compute 0™ ix (vd) through 



Eq. (JlOj) . The first step is to determine the stability cri- 
terion for individual surface traps as a function of their 
shape and orientation. The second step is to develop 
an adequate statistical description of the distributions of 
these surface traps as functions of the shape and orien- 
tation parameters identified in the first step. 



A. Stability of a surface trap 

For a surface trap of specified geometry, represented 
by its base triangle, what is the angle to which the pile 
can be tilted until the trap can no longer stably support a 
sphere and the sphere would roll out? In order to answer 
this question, we first need to quantitatively describe the 
geometry of the surface trap with respect to the surface 
of the pile. This is done in the Appendix, where the yaw 
(f>, roll -0 and tilt 9 of a base triangle are defined (See 
Fig. 0.) 

The determination of stability criteria as a function of 
shape, yaw and roll is a straightforward but tedious job. 
We have used a Mathematica notebook to compute the 
stability diagram for equilateral traps as a function of 
normalized average edge length a = (h + h + h)/(3d), 
yaw 0, and roll ip, where d is the diameter of the spheres. 

The dependence of the maximum stability angle 9 ma x 
on a for traps with ip — is plotted in Fig. ga. It is clear 
that this parameter greatly influences the stability of the 
pile. In order to estimate the value of a for a random 
packing of spheres with a given packing fraction let 
us consider the tetrahedra in a Delaunay tessellation of 
the packing. Provided that the number of tetrahedra per 
sphere do not change for the packings of interest, the av- 
erage volume of the tetrahedra will vary as Vtet ~ 3> _1 . 
Spheres on the surface layer will settle into the minima of 
their traps, thus the tetrahedra they from together with 
their three supporting spheres always have three edges 
whose lengths are equal to the diameter d. Thus, the av- 
erage edge length of the faces that form the base triangle 
is expected to vary as l av ~ $ -1 / 2 . Since all edge lengths 
are equal to d for the densest packing with $ cp = 0.74, 
the estimate for the average normalized edge length of 
base triangles is 



0J4 



-1/2 



(11) 



For the monomer pile with $ = 0.58, this gives a = 1.13, 
in agreement with direct measurements done on the pile 



(see Sec. [VBj.) 

Fig. ||b shows 9 max against and -0 for surface traps 
with a = 1.13. (For certain values of yaw and roll, there 
is no tilt angle for which the traps are stable, and there- 
fore 9 max is undefined.) 

The analysis in Ref. Q is more limited in the types of 
surface traps it considers, as it only looks at traps with 
a = 1 and = 0. As seen in Fig. roll and edge length 
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have great potential impact on the stability of a surface 
trap. Although the Mathematica notebook can deter- 
mine the stability diagram for the most general case, we 
will restrict our analysis to equilateral traps in order to 
keep the subsequent analysis tractable. 



V m (0) 




FIG. 8. The maximum stability angle 8 m ax of equilateral 
traps as a function of (a) normalized edge length a for roll 
■0 = 0, and yaw <j> — (solid line) and (f> — 30° (dashed line); 
(b) 9 m ax as a function of yaw <f> and roll ip for a = 1.13. 



B. Statistical description of traps on the surface of a 

pile 

Having characterized and obtained stability criteria for 
a given surface trap, the next task is to obtain a statis- 
tical description of their population, through probability 
density functions (PDFs). In this study, we will neglect 
short-range correlations between adjacent traps, e.g., as- 
sociated with the sharing of edges, and assume that they 
are drawn independently from an ensemble described by 
PDFs for the values of edge lengths, yaw, roll, and tilt. 
For the present, we will assume that all traps are equi- 
lateral triangles with edge length ad, with uniformly dis- 
tributed yaw angles and Gaussian roll and tilt angle dis- 
tributions: 

'Pahape(h,h,h) = Hh ~ ad)S(l 2 - ad)S(l 3 - ad), (12) 
' (13) 



yawO) = o - ' -TT <(j)<n, 



VroaW 



1 



(14) 



(15) 



In the above, {k} correspond to the edge lengths; we will 
neglect the variability in the shapes of the traps and fo- 
cus on equilateral traps of uniform size in order to study 
the effect of orientational disorder. If desired, the subse- 
quent analysis can be generalized to study the impact of 
disorder in the shapes of the traps as well. 

The orientational PDFs are motivated by assuming 
that the pile surface was created with no initial tilt, and 
rotationally isotropic in the plane of the surface, and that 
little or no rearrangement took place in the surface traps 
during the subsequent tilting of the pile. This would re- 
sult in a uniform PDF of yaw, and nearly identical PDFs 
for roll and tilt (er^ « a$) Jl3| |. 

Fig. |^ depicts how 8 C for a monomer pile changes as a 
function of change in (a) the trap size parameter a, and 
(b) the standard deviations of roll and tilt distributions, 
both individually and jointly. From these plots, it be- 
comes clear that we need additional information about 
the grain-scale roughness of the surface in order to quan- 
titatively predict 9 C for the monomer-dimer piles. 




1.2 



1.4 

a 



1.6 



25 



tilt only- 



roll and tilt 




3 6 9 12 15 

Standard Deviation HdegreesL 

FIG. 9. The dependence of C for a monomer pile as sur- 
face properties are changed from a randomly oriented, flat, 
close-packed surface with no roll, (a) An increase in the nor- 
malized edge length a for equilateral traps stabilizes the traps 
and increases 8 C - (b) Individual effects of including a Gaus- 
sian roll distribution (dotted line, destabilizing), tilt distri- 
bution (dashed line, stabilizing) and the combined effect of a 
simultaneous roll and tilt distribution with the same standard 
deviation (solid line, either stabilizing or destabilizing). 



8 



3 6 

o 
u 




25 



20 
ti 15 

3 
O 

u 10 



160 - 80 



80 160 



f HdegreesL 



3 10 




- 160 


- 80 


80 


160 




f 


HdegreesL 





30 
25 
20 
15 
10 




□DL 



40 - 20 



20 40 



40 



y HdegreesL 



20 20 
y HdegreesL 



40 




35 
30 
25 
20 
15 
10 
5 




D 



D 



- 20 20 40 60 - 10 10 30 50 70 90 

q HdegreesL q HdegreesL 

FIG. 10. Histograms of yaw, roll and tilt distributions for the piles shown in Fig. H. Left: monomer pile, Right: dimer pile. 



In order to test whether real surfaces of piles exhibit 
the assumed behavior, and to obtain representative val- 
ues for the average trap size and the width of yaw, roll 
and tilt distributions, we imaged a portion of the surface 
of a monomer and a dimer pile (see Sec. |l|). The shape 
and orientation of surface traps were identified as follows: 
After locating the centers of the particles on the surface 
layer by stereographic imaging, we computed the average 
surface of the plane by a least square fitting of the cen- 
ters of mass to a plane. We then performed a Delaunay 
triangulation of the particles projected on to this plane 
in order to identify all base triangles associated with po- 
tential surface traps. We then measured the yaw, roll 
and tilt of all the base triangles and created histograms. 
We observed a uniform yaw distribution within a charac- 
teristic sampling error, justifying the use of Eq.(|l3|). We 
also determined the standard deviations {ex,/,, erg} for the 
roll and tilt histograms. The histograms are shown in 
Figs, [l^ for the monomer and dimer pile. 

The comparison between the monomer and dimer piles 
revealed a moderate increase of er^ and erg from about 



15° to 20°, indicating a roughening of the surface along 
with the originally observed reduction in packing frac- 
tion. The average edge length increased from 1.13 to 1.18, 
in agreement with Eq. (pi]). Due to the modest changes 
in these parameters, we have used the trap characteris- 
tics obtained from the monomer pile in the computation 
of 8 C for all the mixture piles. Integrating the stability di- 
agram shown in Fig. with the PDFs given in Eqs.(|l2|- 
|l5| ) to obtain the appropriate DOS{9) through a general- 
ized form of Eq.(^), we finally com put e 9™™{vd) through 
Eq.([|). The result is shown in Fig. [ll]as a solid line, and 
agrees well with experiment for < vj, < 0.6, particularly 
considering that all the parameters have been provided 
by independent measurement. No adjustable parameters 
remain in the model, suggesting that the assumption of 
monomer failure at the surface is valid in this range. The 
disagreement at larger Vd is not surprising, given that 
the pile is comprised almost entirely of dimers and the 
assumption of monomer failure in the theory is expected 
to break down in this limit, resulting in an over-estimate 
of the stability of the pile. 








0.2 0.4 0.6 0.8 1. 
n d 

FIG. 11. Measured and calculated values of tan# c as a 
function of dimer weight fraction Ud- Circles: Experiment. 
Solid Line: Computation based on the statistical description 
of surface traps given by Eqs.(|l^)-(|lE|). Dashed Line: Com- 
parison to the earlier result based on a randomly oriented, 
flat, close-packed surface, reproduced from Fig. m. 



VI. CONCLUSION 

By introducing dimer grains in a sandpile comprised 
of rough spherical monomer grains, we have shown that 
the critical angle of the sandpile can be nearly doubled 
in the limit of high dimer content. Qualitatively, this re- 
sult is not so surprising, given the significant number of 
previous measurements of the angle of repose of mixtures 
of cylindrical or spheri-cylindrical grains with spherical 
grains that have also shown an increase. However, the 
use of dimers, rather than other elongated objects, per- 
mits the random surface of the pile to be described as a 
collection of spheres that form triangular surface traps. 
These triangles have a distribution of edge lengths and 
yaw, roll, and tilt angles that can be directly obtained 
through stereo imaging and can be included in a theory. 
By comparing the macroscopic average property of the 
pile, the critical angle, with the grain-scale structure on 
the pile's surface obtained through imaging, we are able 
to show which aspects of the surface structure are impor- 
tant for determining the value of the critical angle. 

For instance, the treatment of the critical angle of a 
random pile that considers only the mean angle of sta- 
bility for a grain on a close packed surface, averaged over 
the yaw angle, may give a value close to the measured 9 C , 
but is this just a fortuitous agreement? Our results show 
that, by averaging over realistic distributions of yaw and 
tilt, the more realistic median critical angle drops below 
the observed 8 C . However, because the pile is random, 
the intergrain separation has a distribution itself, and 
the average edge length of a triangular surface trap is 
slightly greater than the grain diameter. This increase 
in the edge length of the surface trap, as compared to a 
perfectly close packed surface, increases the stability of a 
sphere in the trap. Indeed, we believe that it is the com- 
bination of the destabilizing influence of the roll distri- 
bution, along with the stabilizing influence of the larger 
edge length that gives the random pile of rough spherical 



grains its rather well-established value of 9 C ss 23°. 

These results for mixed monomer- dimer sandpiles shed 
some light on the observed initial increase in the criti- 
cal angle of wet sandpiles, independent of container size 
and liquid surface tension, when liquid content is below 
a threshold value. @||. The initial linear increase in 
tan 9 C with Vd may provide a plausible mechanism, in 
which the fraction of strongly wetted intergrain contacts 
increases gradually until all intergrain contacts are nearly 
uniformly wetted. Provided that the formed bonds are 
strong enough and relatively dilute, the wet pile can be 
expected to respond similarly to a pile with a small frac- 
tion of dimers. One would expect that clusters with 
grains having more than one cohesive contact with neigh- 
boring grains would form in the wet sandpile as the 
threshold volume fraction is approached, and, above the 
threshold volume fraction, the picture of an average co- 
hesive force holding grains together everywhere in the 
pile would become tenable. The data in Refs. and Q 
suggest an equivalent Ud of about 0.12 at the threshold 
liquid volume fraction. 

It may be possible to extend the presented work to 
systems involving trimers and higher order clusters of 
grains G3]. However, such clusters can have many dif- 
ferent shapes, and to simplify the theoretical treatment, 
it may be necessary to restrict allowed shapes to close- 
packed or linear structures. Along a different direction, 
reducing the grain-grain friction coefficient /i will allow 
sliding failure modes and thereby lower the critical angle 
of the sandpile. Densification of the pile through tapping 
might also change the angles of stability. Finally, de- 
veloping a theoretical understanding of the reduction of 
the grain packing fraction with increasing dimer content 
would help shed light on how strong intergrain attractive 
forces can alter the bulk structure of a random pile. 

We thank P. Chaikin, Z. Cheng, G. Grest and P. Schif- 
fer for stimulating discussions and suggestions. AJL was 
supported in part by the National Science Foundation 
under Award DMR-9870785. 



APPENDIX: PARAMETERIZATION OF A 
SURFACE TRAP - YAW, ROLL AND TILT 

In this Appendix, we define the parameters that de- 
scribe the shape and orientation of a base triangle that 
connects the centers of mass of the three supporting 
spheres that form a surface trap. 

The geometry is shown in Fig. [l^. The coordinate sys- 
tem is fixed such that gravity is in the —z direction, and 
the pile, whose mean surface is initially in the xy— plane, 
is "tilted" by rotating it around the a;— axis. Since overall 
translation of the triangle has no effect on trap stability, 
the vertex across the shortest edge (with length l\) has 
been arbitrarily placed on the 2— axis for ease of illus- 
tration. The base triangle can be fully specified by the 
positions of the two remaining vertices relative to the first 
one; this leaves six parameters to be determined. 
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and therefore they will have independent probability dis- 
tributions. Splitting the parameters that describe these 
two attributes avoids dealing with joint probability dis- 
tributions across these two classes of parameters. 

(ii) Tilting the pile does not change the yaw and roll of 
a surface trap. Thus, a "stability interval" [0 m i n ,9 max ] 
can be defined for a surface trap of given yaw and roll, 
corresponding to all the values of tilt 9 for which the trap 
can stably support a surface particle. 

□ □ 



FIG. 12. The geometry of a surface trap, characterized by 
its yaw <f>, roll tp and tilt 6, as well as the edge lengths {h} of 
its base triangle. 

Given these parameters, the original base triangle can 
be reconstructed (modulo translations) as follows: Place 
a triangle with given edge lengths in the xy— plane such 
that the shortest edge is parallel to the x— axis and 
"downhill" from the vertex across it, i.e., the y— ordinate 
of the vertex is larger. Then, rotate the triangle around 
the z— axis by <f>, y— axis by tp and finally, x— axis by 0. 
With the proper labeling of vertices as described above, 
every triangle is uniquely identified except for degener- 
ate cases (isosceles and equilateral triangles), in which 
case the stability criteria are identical and the particular 
choice of angles is immaterial. 

This parameterization has two main advantages, both 
of which facilitate statistical averaging over many traps, 



performed in Sec. VB: 

(i) The shapes and orientations of surface traps are 
very likely to be statistically independent of each other, 
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